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In 2D electronic spectroscopy stndies, long-lived quantnm beats have recently been 
observed in photosynthetic systems, and several theoretical studies have suggested 
that the beats are produced by quantum mechanically mixed electronic and vibra¬ 
tional states. Concerning the electronic-vibrational quantum mixtures, the impact 
of protein-induced fluctuations was examined by calculating the 2D electronic spec¬ 
tra of a weakly coupled dimer with the Franck-Condon active vibrational modes 
in the resonant condition [Fujihashi et ah, J. Chem. Phys. 142, 212403 (2015), 
arXiv:1505.05281,]. This analysis demonstrated that quantum mixtures of the vi- 
bronic resonance are rather robust under the influence of the fluctuations at cryo¬ 
genic temperatures, whereas the mixtures are eradicated by the fluctuations at phys¬ 
iological temperatures. However, this conclusion cannot be generalized because the 
magnitude of the coupling inducing the quantum mixtures is proportional to the inter¬ 
pigment electronic coupling. In this study, we explore the impact of the fluctuations 
on electronic-vibrational quantum mixtures in a strongly coupled dimer with an off- 
resonant vibrational mode. Toward this end, we calculate energy transfer dynamics 
and 2D electronic spectra of a model dimer that corresponds to the most strongly 
coupled bacteriochlorophyll molecules in the Fenna-Matthews-Olson complex in a 
numerically accurate manner. The quantum mixtures are found to be robust under 
the exposure of protein-induced fluctuations at cryogenic temperatures, irrespective 
of the resonance. At 300 K, however, the quantum mixing is disturbed more strongly 
by the fluctuations, and therefore, the beats in the 2D spectra become obscure even 
in a strongly coupled dimer with a resonant vibrational mode. Further, the over¬ 
all behaviors of the energy transfer dynamics are demonstrated to be dominated by 
the environment and coupling between the 0 — 0 vibronic transitions as long as the 
Huang-Rhys factor of the vibrational mode is small. The electronic-vibrational quan¬ 
tum mixtures do not necessarily play a signihcant role in electronic energy transfer 
dynamics despite contributing to the enhancement of long-lived quantum beating in 
the 2D spectra. 
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I. INTRODUCTION 


The development of ultrafast nonlinear spectroscopic techniques has provided an incisive 
tool to probe the quantum dynamics of condensed phase systems.^ Investigating photo¬ 
synthetic pigment-protein complexes using two-dimensional (2D) electronic spectroscopy 
revealed long-lived beating phenomena in the 2D spectrspHIIl Engel et alP revealed the 
presence of quantum beats persisting for at least 660 fs in the 2D electronic spectra of the 
Fenna-Matthews-Olson (FMO) pigment-protein complex from green sulfur bacteria at cryo¬ 
genic temperatures. Panitchayangkoon et ah® subsequently demonstrated that quantum 
beats in the FMO complex persist for at least 1.5 ps and 300 fs at cryogenic and physio¬ 
logical temperatures, respectively. Recently, nuclear vibrational effects have been theorized 
to explain the spectral beating behavioi^^^®^, which persist for significantly longer than the 
predicted electronic dephasing times.^^^^ Christensson et al.^ proposed that the resonance 
between electronic states and the Franck-Condon active vibrational states serves to create 
vibronic excitons, i.e., quantum mechanically mixed electronic and vibrational states. Such 
states have vibrational characteristics and enhanced transition dipole moments owing to 
intensity borrowing from strong electronic transitions, despite the small Huang-Rhys factors 
of bacteriochlorophyll (BChl).l23S6| Therefore, vibronic resonance was suggested to produce 
oscillations in the 2D signal that persist for several picoseconds. In this regard, Tiwari, 
Peters, and Jona^ showed that vibronic resonance also enhances the vibrational coherence 
in the electronic ground state and demonstrated that this effect could explain the long-lived 
oscillations in the 2D spectra. 

Indeed, the concept of vibronic resonance plausibly explains the spectroscopically ob¬ 
served long-lived spectral beating. However, the dynamic interplay between the electronic- 
vibrational quantum mixture and electronic energy fluctuations induced by the environmen¬ 
tal dynamics naturally raises a questionj ^^ * ^^ ! To address this question, Fujihashi, Fleming, 
and IshizakP^ considered a weakly coupled model dimei^ inspired by BChls 3 and 4 in the 
FMO complex of Chlorobaculum tepidum {C. tepidurn^^, in which the gap between the two 
electronic energy eigenstates resonates with one of the strongest vibrational modes with fre¬ 
quency 180 cm~^ in the absence of the surrounding environment. They numerically explored 
the impacts of fluctuations on the quantum mixture of the vibronic resonance by analyzing 
the beating amplitudes of the ground state bleaching contribution in calculated nonrephas- 
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ing 2D spectra. They demonstrated that the vibronic resonance is rather robust under 
the influence of environmentally induced fluctuations at cryogenic temperatures, resulting 
in long-lasting beats of vibrational origin in the 2D spectra. However, larger-amplitude 
fluctuations disturbed the quantum mixture of the vibronic resonance more strongly at 
physiological temperatures, causing the beats in the 2D spectra to disappear. Further, they 
suggested that vibronic resonance does not play a significant role in the electronic energy 
transfer (EET) dynamics, at least in the weakly coupled dimer, despite contributing to the 
enhancement of the long-lived quantum beating in the 2D electronic spectra. 

However, this conclusion cannot be generalized. Here, it should be noted that the mag¬ 
nitudes of the couplings inducing electronic-vibrational quantum mixtures are dependent 
on inter-pigment electronic couplings as well as the Huang-Rhys factors of the vibrational 
modes. Indeed, there exist several pairs of BChls constituting strongly coupled dimers in the 
FMO complex. Even if the gaps between the electronic energy eigenstates of such strongly 
coupled dimers do not resonate with relatively strong vibrational modes in BChl molecules, 
it may be possible to induce relevant electronic-vibrational quantum mixtures that affect the 
time evolution of the 2D spectra and/or corresponding EET dynamics. This study mainly 
aims to explore the influences of quantum mechanically mixed electronic and vibrational 
states in a strongly coupled model dimer under the exposure of protein-induced fluctua¬ 
tions. Toward this end, we study the 2D electronic spectra and EET dynamics in a model 
dimer mimicking BChls 1 and 2 of the FMO complex, which is the most strongly coupled 
pair. 


II. THEORY 

We consider a coupled hetero-dimer, the simplest EET system. To describe the EET 
dynamics and optical responses, we restrict the electronic spectra of the m-th pigment in 
a pigment-protein complex (PPG) to the ground state, and the first excited state, 

\^me)) although higher excited states are occasionally of consequence in nonlinear spectro¬ 
scopic signals.®^ No experimental evidence has been reported for nonadiabatic transitions 
and radiative/nonradiative decays between \ipme) and \(pmg) in light-harvesting complexes on 
picosecond timescales. Thus, we organize the product states in the order of their elementary 
excitation numbers. The overall ground state with zero excitation reads \g) = \(pig)\(p 2 g)■ 
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The presence of a single excitation at pigment 1 is expressed as |1) = \^ie)\^ 2 g) and the 
other, as |2) = \(pig)\(p 2 e) ■ As the intensity of sunlight is weak, the single-excitation manifold 
is of primary importance under physiological conditions. However, nonlinear spectroscopic 
techniques such as 2D electronic spectroscopy can populate some higher excitation mani¬ 
folds, e.g., the double-excitation manifold comprising |12) = \Lpie)\(p 2 e)■ Therefore, the PPG 
Hamiltonian can be expressed as 

hfppc ~ Hgxc T hfgnv T hfgxc —env? (2.1) 

The hrst part Hexc is fhe Hamiltonian describing the electronic transition expressed as 

iPexc = Eg\g){g\ 

2 2 
+ J2iEg + hnm)\m){m\ + J2Y1 hJmn\m){n\ 

m=l m=l n^m 

+ {^Eg + hQi + /ir22)|12)(12|, (2-2) 

where hQm is the Franck-Con don transition energy of the mth pigment, and hJmn denotes 
the electronic coupling between the mth and the nth pigments. In the derivations above, 
the energy of the overall electronic ground state. Eg, has been introduced; however, we 
set Eg = 0 without loss of generality. The second term Henv is the Hamiltonian of the 
environmental and vibrational DOFs. The last term hfgxc-env describes the coupling of the 
environmental and vibrational degrees of freedom (DOFs) to the electronic excitations as 

2 

Hexc—env ni\rn){m\ + {ui + U2)\12){12\. (2.3) 

m=l 

The operator Um. is termed the collective energy gap coordinate to describe the fluctuations 
in the electronic energy and the dissipation of the reorganization energy. In this study, we 
assumed that the fluctuations in the electronic energies of different pigments are not corre¬ 
lated. It should be emphasized that the coordinate, Um, includes information associated with 
not only the electronic excited state but also the electronic ground state, as demonstrated 
in Ref.^. The electronic coupling between the pigments may also depend on the environ¬ 
mental and nuclear vibrational DOFs. However, in this study, we assume that the nuclear 
dependence of hJmn is vanishingly small, and we employ the Condon-like approximation as 
usual. 
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We assume that the environmentally induced fluctuations can be described as Gaus¬ 
sian processes and that the relevant nuclear dynamics are harmonic vibrations. Under this 
assumption, the dynamic properties of the environmental and intramolecular vibrational mo¬ 
tions are characterized by the two-body correlation functions of Um{t) = 
as 


(Mm(f)Mm(O)) = D^{t) - 


(2.4) 


where the canonical average has been introduced, (...) = Tr(... e ''^'^®"'')/Tr e with (3 


being the inverse temperature, l/ksT. In Eq. (2.4), Dmif) is the symmetrized correlation 
function of itmit), defined as Dmit) = (1/2) and $m(f) is the response 

function defined by <hm(^) = (V^) {[umit), Um{0)]-)■ The quantum fluctuation-dissipation 
relatiorP^ allows one to express the symmetrized correlation and response functions as 


D„(t) = - 

71 

2 

^m(t) = — / duJm{u}) sinojt, (2.6) 

^ Jo 

where Jm{^) is termed the spectral density and is defined by the imaginary part of the 
Fourier-Laplace transform of the response function. The absolute magnitude of the spectral 
density is related to the total reorganization energy, hXm,env + h\m,vih = duJm{u})/{7lu), 

where hXm,env and hXm,vih denote the environmental and vibrational reorganization ener¬ 
gies of the mth pigment, respectively. The spectral density can be separated into two- 
components: a continuous part corresponding to the environmental contribution and discrete 
peaks corresponding to the vibrational modes, namely, P^ * ^^ * ^^ * ^^ * ^ 


' r / \ 1 

dUjJm\, 0 O) coth —— cos ut, 


(2.5) 




J n 


r(^) T Jm,vib(^) ■ 


(2.7) 


To focus on the timescale of the environmental dynamics affecting the electronic transition 
energies, we model the environmental component as the Drude-Lorentz spectral density as 


Tm,env(^) 2hXg 


'Xenv^ 

env 


' 9 I 

+ 7, 


( 2 . 8 ) 


where corresponds to the timescale of the environment-induced fluctuations. For sim¬ 
plicity, we consider a single intramolecular vibration on each pigments with frequency cUvib 
and Huang-Rhys factor S. We describe the spectral density for the vibrational mode using 
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the Brownian oscillator model with the vibrational relaxation rate 7vib, such that 


'-^m,vib(^) ^-^vib 




vib 


7vib 


7vib 


(2.9) 


Wvib L(^-^vib)^ + 7vib (^ + ‘^vib)^ + 7vib. 

where oivib = (i^vib “ 7vib)^'^^ ^-^vib = li^vihS. In the slow vibrational relaxation limit 
(7vib —t 0), Eq. ( |2.9[ ) is reduced to Jm,vih{^) = '7ihSu‘^[6{u — Wvib) — 5{u + Wvib)], which has 

been employed in the literature .^^SES] 

An adequate description of the EET dynamics is given by the reduced density operator, 
Pexc(^), that is, the partial trace of the total PPG density operator, pppc(^), over the envi¬ 
ronmental and nuclear DOFs: Pexcit) = TrenvPppc(^)- For this reduction, we assume that 
the total PPG system at the initial time {t = 0) is in the factorized product state of the form 
Pppc ~ Pexc®e“^^™''/Tr In other words, Pexc(^) = Trenv[G(^)Pppc]) where G{t) is the 

retarded propagator of the total PPG system in the Liouville space. This factorized initial 
condition is generally considered to be unphysical in the literature on open quantum systems 
as it neglects the inherent correlation between a system of interest and its environment 
However, in electronic excitation processes, this initial condition is inconsequential because 
it corresponds to the electronic ground state or electronic excited state generated in accor¬ 
dance with the vertical Franck-Gondon transition.^^^® The Gaussian property of Um enables 
one to derive a formally exact equation of motion that can describe the EET dynamics un¬ 
der the influence of the environmentally induced fluctuation and dissipation. Therefore, the 
EET dynamics under the influence of the environmental/nuclear dynamics and correspond¬ 
ing 2D electronic spectra can be described in a numerically accurate manner. The technical 
details of the computations are given in 


III. RESULTS AND DISCUSSION 

In this section, we present and discuss the numerical results to explore the influences of 
quantum mechanically mixed electronic and vibrational states in a strongly coupled dimer 
under the exposure of protein-induced fluctuations. Toward this end, we address a model 
dimer mimicking BGhls 1 and 2 in the FMO complex of C. tepidun^^^^, this is the most 
strongly coupled pair and exhibits wave-like EET persisting for several hundred femtoseconds 
even in the absence of vibronic contributions.^ The Franck-Gondon transition energies of 
BGhls 1 and 2 are set to Di = 12410 cm“^ and D 2 = Di -|- 120 cm“^, and their electronic 
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coupling is set to J 12 = —87.7cm ^ according to Ref.SS!. The transition dipole moment 
directions are set to = 66.6° based on the PDB hie 3ENI.®I 

In the descriptions below, \xTv) denotes the n-th vibrational level of the a-th electronic 
state in the m-th pigment. For simplicity, we consider a single vibrational mode with 
frequency cjvib of BChl 1 only, and we focus on the quantum mixture between the vibronic 
transition, |9?ig)|Xgo) 17^16)1X61)) of BChl 1 and the electronic transition, \(p 2 g) -H- \^ 2 e)) 

of BChl 2. The mixture is induced by the coupling ^i2(XeilXgo) ~ — J 12 exp(—S'/2)\/5, 
where S denotes the Huang-Rhys factor of the mode. The vibrational modes of BChl in 
solutions or protein environments have been investigated,^^^® and the strongest vibrational 
mode has been found at approximately 180 cm“^. However, this vibrational mode does 
not resonate with the gap between the electronic energy eigenstates in the model dimer 
under investigation, [(f22 — + 47 ^ 3 ]^'^^ = 213cm“^. Therefore, we investigate two cases 

concerning the vibrational frequency: the off-resonance case, cUvib = 180 cm“^, and the 
resonance case, cUvib = 213 cm“^. 

For the numerical calculations of the 2D spectra and FET dynamics, the environmental 
reorganization energy and reorganization time constant are set to Aenv = 35 cm“^ and 7 ^^^ = 
100 fs, and the vibrational relaxation rate is set to y^jb = 2 ps. The Huang-Rhys factors of 
BChl are typically small,^ and the factor associated with the vibrational mode in this work 
is set to S' = 0.025,^ which is within the range of the experimentally measured values.®^ 

A. Beating in 2D electronic spectra 

In this subsection, we explore the contributions of the vibrational modes of the strongly 
coupled dimer in 2D electronic spectra. In Fig. [T| we consider 2D spectra in the ab¬ 
sence of any vibrational contributions as a reference. Figure Eti shows the diagonal cut 
of the calculated nonrephasing 2D spectra as a function of the waiting time at 77 K, 
and Fig. shows a plot of the time evolution of the amplitudes of the diagonal peaks 
(cui = uj^ = 12343cm“^, 12582cm“^). The contour plots in Fig. apparently indicate 
the monotonous decay of the upper diagonal peak induced by the exciton relaxation as 
well as the monotonic increase of the lower diagonal peak, which is the inherent nature of 
nonrephasing signals irrespective of the EFT dynamics. Nevertheless, Fig. reveals the 
presence of spectral beating persisting for 600 fs. This quantum beating originates from the 




strong electronic coupling between BChls 1 and 2 and is thus of electronic origin. Indeed, 
this electronic coherence and its lifetime are consistent with the previous theoretical anal¬ 
ysis of EET dynamics in the FMO complex.^ Figure shows the time evolutions of the 
calculated nonrephasing 2D spectra at 77 K with the vibrational mode of BChl 1 for (a) 
the off-resonance case, Wvib = 180cm“^, and (b) the resonance case, cUvib = 213cm“^. The 
spectral beats of the 2D spectra are enhanced by the vibrational contributions and persist 
for at least 2 ps. This observation demonstrates that even in the off-resonance case, the 
electronic-vibrational quantum mixture causes distinct beats in the 2D spectra of strongly 
coupled dimers at cryogenic temperatures. 

However, the variance of the amplitude of the environmentally induced fluctuations in the 
electronic energy is estimated via Eq. (2.5) as Dm,env(0) — 2 /cbT /nu, where 

the high-temperature approximation, coth(/3ha;/2) ~ 2kBT/hw, has been employed. The 
fluctuations at higher temperatures are larger in magnitude, and therefore, the electronic- 
vibrational quantum mixture is supposedly fragile at physiological temperatures. Figure]^ 
shows the time evolution of the diagonal cut of the calculated nonrephasing 2D spectra 
at 300 K in the absence of the vibrational contribution, and Figs. i) and 1^ show the 
2D spectra with the vibrational mode in the off-resonance case, Wvib = 180 cm“^, and the 
resonance case, cUvib = 213 cm“^, respectively. In both the off-resonance and resonance 
cases, the amplitudes of the beating decrease signihcantly, although Figs. it and show 
spectral beating arising from the electronic-vibrational mixture, albeit only slightly. This 
indicates that the environmentally induced fluctuations at physiological temperatures almost 
destroy the quantum mixing even in a strongly coupled dimer. In this study, inhomogeneous 
broadening in the 2D spectra that is caused by the static disorder of the electronic transition 
energies has not been considered. Under the influence of inhomogeneous broadening, as seen 
in the experimental data, the small amplitudes of the beats observed in Fig.j^are presumably 
disturbed and are difficult to observe. 


B. Population dynamics 

We explore the influence of intramolecular vibration on the spatial dynamics of electronic 
excitation. Figure shows the time evolution of the energy donor (BChl 2) population 
affected by the intramolecular vibration at temperatures of (a) 77 K and (b) 300 K. The 
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red and green lines indicate the dynamics in the presence of vibrations with freqnency 
cjvib = 180 cm“^ and cUvib = 213 cm“^, respectively. The bine lines indicate the dynamics 
withont the vibrations as a reference. The EET dynamics in Fig. |^show wave-like motion 
induced by the quantum superposition between the 0—0 vibronic transition of BChl 1 and the 
electronic transition of BChl 2 because of the strong electronic coupling. In the case without 
a vibrational mode, the electronic coherence persists for 700 fs and 300 fs at 77 K and 300 
K, respectively. The coupling constant inducing the electronic-vibrational quantum mixture 
'^i2(XeilXgo) proportional to the electronic coupling. Consequently, the EET dynamics 
in the case of strong electronic coupling shows long-lasting (up to at least 2 ps) wave-like 
motion induced by the quantum mechanically mixed electronic and vibrational states in 
contrast to that in the case of the small electronic coupling investigated in Ref.^. The 
vibrational mode slightly contributes to the acceleration of the EET because the electronic- 
vibrational quantum mixture adds a new EET channel originating from the 0 — 1 vibronic 
transition of BChl 1 and the electronic transition of BChl 2. As the vibrational frequency 
approaches the resonance condition, the contribution of the vibrational mode is increased, 
as shown in Fig. However, the overall behaviors of EET dynamics are dominated by 
the environment and the coupling between the 0-0 vibronic transition of BChl 1 and the 
electronic transition of BChl 2. The inter-site coupling constant inducing the the quantum 
mixture, Ji 2 {x\i\x])o) = —T 12 exp(—S'/2)A/iS = 13.7cm“^, is much smaller than the coupling 
constant between the 0-0 vibronic transition of BChl 1 and the electronic transition of 
BChl 2, Ji 2 {x\o\x\o) — <^12 exp(—5/2) = —86.6cm“^, which gives rise to the fast channel, as 
long as the Huang-Rhys factor is small. Consequently, the electronic-vibrational quantum 
mixture does not play a signihcant role in enhancing the EET dynamics under the influence 
of the environmentally induced fluctuations even in the most strongly coupled BChls in the 
FMO complex, despite contributing to the enhancement of the long-lived beating in the 2D 
electronic spectra. 


IV. CONCLUDING REMARKS 

In this study, we have explored the impact of environmentally induced fluctuations on 
the quantum mechanically mixed electronic and vibrational states of pigments in a strongly 
coupled dimer. Toward this end, we investigated BChls 1 and 2, the most strongly coupled 
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dimer in the FMO complex, and performed accurate numerical calculations of the EET 
dynamics and corresponding 2D electronic spectra. 

In the strongly coupled dimer at cryogenic temperatures, the electronic-vibrational quan¬ 
tum mixing is rather robust against the environmentally induced fluctuations irrespective of 
the resonance, resulting in long-lasting spectral beating of vibrational origin in the 2D elec¬ 
tronic spectra. However, the amplitudes of the environmentally induced fluctuations in the 
electronic excitation energy increase with the temperature. At physiological temperatures, 
the electronic-vibrational quantum mixture is disturbed more strongly, and therefore, the 
beating in 2D spectra becomes obscure despite the strong electronic coupling. 

Furthermore, we have investigated the effect of the underdamped intramolecular vibra¬ 
tional motion on the EET dynamics using the same equation of motion and parameters as 
in the calculations of the 2D electronic spectra. The electronic-vibrational quantum mix¬ 
ture leads to an oscillation lasting up to several picoseconds at cryogenic and physiological 
temperatures. However, the overall behavior of the EET dynamics is dominated by the envi¬ 
ronment and the coupling between the 0 — 0 vibronic transition of BChl 1 and the electronic 
transition of BChl 2 as long as the Huang-Rhys factor is small. The electronic-vibrational 
quantum mixture contributes only slightly to the acceleration of the electronic energy trans¬ 
fer, although the 2D electronic spectra clearly exhibit long-lived beating enhanced by the 
quantum mixture. This suggests that information about the oscillatory behavior extracted 
from the 2D spectra should be interpreted carefully. 
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FIG. 1. (a) Time evolution of the diagonal cut of the nonrephasing 2D spectra at 77 K of a coupled- 
dimer without intramolecular vibration, (b) Time evolutions of the amplitudes of the lower diagonal 
peak, (wi = UJ 3 = 12,343cm“^), and the upper diagonal peak, (wi = ws = 12,582cm“^), in the 
spectra of Fig. The environmental reorganization energy and the reorganization time constant 
are set to Aenv = 35cm“^ and 7 ^^ = 100 fs, and the Huang-Rhys factor and the vibrational 
relaxation time constant are S = 0.025 and = 2ps. The normalization of the plots is such that 
the maximum value of the nonrephasing 2D spectra of Fig. at t 2 = Ofs is unity, and equally 
spaced contour levels (0, ±0.01, ±0.02,...,) are drawn. To clarify the relationship between the 
beating amplitude magnitudes in Fig. [3>, the peak amplitudes are shifted. 
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FIG. 2. Time evolution of the diagonal cut of the nonrephasing 2D spectra at 77 K of a coupled- 
dimer with intramolecular vibration of frequency (a) Wvib = 180 cm“^ and (b) Wvib = 213 cm“^. 
The normalization of the plots is such that the maximum value of the nonrephasing 2D spectra in 
Fig.[^ at t 2 = Ofs is unity, and equally spaced contour levels (0, ±0.01, ±0.02,...,) are drawn. 
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FIG. 3. Time evolution of the diagonal cut of the nonrephasing 2D spectra of (a) a coupled-dimer 
without intramolecular vibration, (b) a coupled-dimer with intramolecular vibration of frequency 
(^vib = 180cm“^, and (c) a coupled-dimer with vibration of frequency cjvib = 213cm“^ at 300 K. 
The normalization of the plots is such that the maximum value of the nonrephasing 2D spectra in 
Fig.^ at t 2 = Ofs is unity, and equally spaced contour levels (0, ±0.01, ±0.02 ,...,) are drawn. 
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FIG. 4. Time evolution of the energy donor population affected by intramolecular vibration of BChl 
1 (acceptor) at temperatures of (a) 77 K and (b) 300 K. In each panel, the red and green lines indi¬ 
cate the dynamics in the presence of vibrations of frequency tUvib = 180 cm“^ and Wvib = 213 cm“^, 
respectively. The blue lines indicate the dynamics without the vibrations. The calculations were 
performed with the parameters employed in Fig. 
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